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For the first time, a next-to-leading BFKL study of the cross section and azimuthal decorrellation 
of Mueller Navelet jets is performed, i.e. including next-to-leading corrections to the Green's 
function as well as next-to-leading corrections to the Mueller Navelet vertices. The obtained re- 
sults for standard observables proposed for studies of Mueller Navelet jets show that both sources 
of corrections are of equal and big importance for final magnitude and final behavior of observ- 
ables, in particular for the LHC kinematics investigated here in detail. The astonishing conclu- 
sion of our analysis is that the observables obtained within the complete next-lo-leading order 
BFKL framework of the present paper are quite similar to the same observables obtained within 
next-to-leading logarithm DGLAP type treatment. The only noticeable difference is the ratio the 
azimuthal angular moments (cos2<p) / (cos <p) which still differs in both treatments. 
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1. Introduction 



The high energy regime of QCD is one the key questions of particle physics. In the semi- 
hard regime of a scattering process in which s>-(, logarithms of the type [a s ln(5/|?|)] n have 
to be resummed, giving the leading logarithmic (LL) Balitsky-Fadin-Kuraev-Lipatov (BFKL) [1] 
Pomeron contribution to the gluon Green's function describing the exchange in the ^-channel. The 
question of testing such effects experimentally then appeared. Based on new experimental facilities, 
characterized by increasing center-of-mass energies and luminosities, various observables have 
been proposed and often tested in inclusive [2], semi-inclusive [3] and exclusive processes [4]. 
The basic idea is to select specific observables in order to minimize standard collinear logarithmic 
effects a la DGLAP [5] with respect to the BFKL one. This aims to choose them in such a way that 
the involved transverse scales would be of similar order of magnitude. 

We will here consider one of the most famous testing ground for BFKL physics: the Mueller 
Navelet jets [6] in hadron-hadron colliders, defined as being separated by a large relative rapidity, 
while having two similar transverse energies. One thus expects an almost back-to-back emission in 
a DGLAP scenario, while the allowed emission of partons between these two jets in a BFKL treat- 
ment leads in principle to a larger cross-section, without azimuthal correlation between them. The 
predicted power like rise of the cross section with increasing energy has been observed at the Teva- 
tron /^-collider [7], but the measurements revealed an even stronger rise than predicted by BFKL 
calculations. Besides, the leading logarithmic approximation [8] overestimates this decorrelation 
by far. Improvements have been obtained by taking into account some corrections of higher order 
like the running of the coupling [9]. Calculation with the full NLL BFKL Green's function [10] 
have been published recently [11]. We here review results of Ref. [12] on the full NLL BFKL 
calculation where also the NLL result for the Mueller Navelet vertices [13] is taken into account. 



2. NLL calculation 



The two hadrons collide at a center of mass energy s producing two very forward jets, whose 
transverse momenta are labeled by Euclidean two dimensional vectors kj j and k/2, while their 
azimuthal angles are noted as <j>j t \ and (pj2- We will denote the rapidities of the jets by yj^ and yj : 2 
which are related to the longitudinal momentum fractions of the jets via xj = ^e yj - Since at the 
LHC the binning in rapidity and in transverse momentum will be quite narrow [14], we consider the 
case of fixed rapidities and transverse momenta. Due to the large longitudinal momentum fractions 
xj.i and x/ 2 of the forward jets, collinear factorization holds and the differential cross section can 
be written as 



da 



— =W dx! / 1 dx 2 / a (x 1 )/ b (x 2 ) d(Ji ; b , (2.1) 

7,1^7,2 tbJo Jo d|k 7i i|d|k 7 , 2 |dj7,id3; 7> 2 



d|ky ; i | d|k^, 2 1 d^y : i dyy, 2 % 

where / a b are the parton distribution functions (PDFs) of a parton a (b) in the according proton 
(which are renormalization scale Hr and factorization scale jif dependent). The resummation of 
logarithmically enhanced contributions are included through kj -factorization: 



dd, 



ab 



d|k A i | d|ky 5 2 1 d^,i dyy,2 



Jd^idfaa y n d 2 k 1 d 2 k 2 y a (-ki,x 1 )G(k 1 ,k2,f)Vb(k2,x 2 ), (2.2) 
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where the BFKL Green's function G depends on s = x\%2S. The jet vertices V a j, were calculated at 
NLL order in Ref. [13]. Combining the PDFs with the jet vertices one writes 



J d0 7il d0 7 , 2 J d 2 kid 2 k 2 <J>(k / , 1 ,^, 1 ,-k 1 )G(k 1 ,k 2 ,f)<I>(k 7i2 ,x /!2 ,k2) 



da 

d|k /; i|d|k 7 , 2 |dj7,idj 7 , 2 

where <J>(k/, 2 ,x 7;2 ,k 2 ) = y<lx 2 /(x 2 )y(k 2 ,x 2 ). (2.3) 



In view of the azimuthal decorrelation we want to investigate, we define the following coefficients: 



= J d0/,i dfe cos - 7j2 - n))Jd 2 1s.i d 2 k 2 <£(k/,i,X/,i , -ki)G(ki,k 2 , j)0(k/, 2 ,xy, 2 ,k 2 ), 



from which one can easily obtain the differential cross section and azimuthal decorrelation as 

^0 and (cos(m<p)) = (cos (m(0y,i - (j)j >2 -7t))) = (2.4) 



d | k /a | d | ky, 2 1 djy, i dyy >2 
The guiding principle is then to rely on the LL-BFKL eigenfunctions 

^, v (k 1 ) = -L(k 2 )' V ^^', (2.5) 
7Tv2 

although they strictly speaking do not diagonalize the NLL BFKL kernel. In the LL approximation, 

f ( s \ m ( m ^ 

tf m = (4-3S m>0 ) J dvC w , v (|k /l i| J ^,i)C^ v (|l^,2|,x /l 2) ( -J , (2.6) 

where 

C m ,v(|k/|,Xy) = J d(t>jd 2 kdxf(x)V{k,x) (k)cos(m0y), (2.7) 

and co(n,v) = a s % Q (|n|, ± + iv) , with Xo(n,f) = 2¥(1) - Y (y+ §) - V (l - y+ f ) and a,. = 
N c a s / 71. The master formulae of the LL calculation ( |2.6| , |2.7[ ) will also be used for the NLL cal- 
culation. The price to pay in order that E ny remains eigenfunction is to accept that the eigenvalue 
become an operator containing a derivative with respect to v. In combination with the impact fac- 
tors the derivative acts on the impact factors and effectively leads to a contribution to the eigenvalue 
which depends on the impact factors. 

At NLL, the jet vertices are intimately dependent on the jet algorithm [13]. We here use the 
cone algorithm, which is expected to the used by the CMS collaboration. At NLL, one should also 

^2 

pay attention to the choice of scale sq. We find the choice of scale sq = v /Jol ^0,2 with sq \ = -j-kj j 

' x j,i ' 

rather natural, since it does not depend on the momenta ki )2 to be integrated out. Besides, the 
dependence with respect to sq of the whole amplitude can be studied, when taking account the fact 
that both the NLL BFKL Green function and the vertex functions are sq dependent. In order to study 
the effect of possible collinear improvement [15, 16], we have, in a separate study, implemented 
for n = the scheme 3 of Ref. [15]. This is only required by the Green function since we could 
show by a numerical study that the jet vertices are free of 7 poles and thus do not call for any 



collinear improvement. In practice, the use of Eqs. (2.6, |2.7|) leads to the possibility to calculate 
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for a limited number of m the coefficients C my as universal grids in v, instead of using a two- 
dimensional grid in k space. We use MSTW 2008 PDFs [17] and a two-loop strong coupling 
with a scale /i« = v/Jltyi | • [k/^l • Although a BFKL treatment does not require in principle any 
asymmetry between the two emmited jets, large logarithms of non-BFKL origin, of Sudakov type, 
could dominate the observable for very symmetric jet configurations. In order to compare our 
analysis with DGLAPNLO approaches [18] obtained through the NLL-DGLAP partonic generator 
Duet [19], for which symmetric configurations lead to instabilities, we here display our results 
for |k/i| = 35GeV, |k/2| = 50GeV, and compare them with the DGLAP NLO prediction (see 
Ref. [12] for symmetric configurations). 

3. Results 

We now present our results for the LHC at the design center of mass energy ^f~s = 14TeV. 
Motivated by a recent CMS study [14] we restrict the rapidities of the Mueller Navelet jets to the 
region 3 < \yj\ < 5, thus limiting the relative rapidity Y between 6 and 10 units. Fig. [I]a and [l|b 
respectively display the cross-section and the azimuthal correlation 1 . This explicitely shows the 
dramatic effect of the NLL vertex corrections, of the same order as the one for the Green function. 
In particular, the decorrelation based on our full NLL analysis is very small, similar to the one based 
on NLO DGLAP. The main source of uncertainties is due to the renormalization scale jj. R and to the 




Figure 1: Differential cross section (a), azimuthal correlation (cos<p) (b) and ratio (cos 2(p) / (cos <p) (c) in 
dependence on Y for |k/,i| = 35GeV, |ky,2| = 50GeV. The errors due to the Monte Carlo integration are 
given as error bands. As dots are shown the results of Ref. [18] obtained with Duet [19]. 



energy scale This is particularly important for the azimuthal correlation, which, when includ- 
ing a collinear improved Green's function, may exceed 1 for small = pLp . The only remaining 
observable for which a noticeable difference can be expected between BFKL and DGLAP type of 
treatment is the ratio (cos 2<p) / (cos <p) , which would deserve a precise experimental analysis. 

In conclusion, contrarily to the general belief, the study of Mueller Navelet jets may not be 
conclusive to exhibit differences between BFKL and DGLAP dynamics. 

Work supported in part by the Polish Grant N202 249235, the grant ANR-06-JCJC-0084 and 
by a PRIN grant (MIUR, Italy). 

'We use the same color coding for both plots, namely blue shows the pure LL result, brown the pure NLL result, 
green the combination of LL vertices with the collinear improved NLL Green's function, red the full NLO vertices with 
the collinear improved NLL Green's function. 
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